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Heavy stable fourth generation Majorana neutrinos contribute to a small fraction of the relic den- 
sity of dark matter (DM) in the Universe. Due to its relatively strong coupling to the standard 
model particles, it can be probed by the current direct and indirect DM detection experiments 
even it is a subdominant component of the halo DM. We show that the current Xenon 100 data 
constrain the mass of the stable Majorana neutrino to be greater than the mass of the top quark. 
The effective spin-independent cross section for the neutrino elastic scattering off nucleon is pre- 
dicted to be ~ 1.5 X 10^'*'* cm^, which is insensitive to the neutrino mass and mixing and can be 
reached by the direct DM detection experiments in the near future. In the same mass region the 
predicted effective spin-dependent cross section for the heavy neutrino scattering off proton is in 
the range of 2 x lO^'*" cm^ ~ 2 x 10^^^ cm^, which is within the reach of the ongoing DM indi- 
rect search experiments. We demonstrate such properties of the heavy neutrino DM in a fourth 
generation model with the stability of the fourth Majorana neutrino protected by an anomaly-free 
U{1) gauge symmetry. 
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Models with chiral fourth generation fermions are among the simplest and well-motivated 
extensions of the standai^d model (SM) and have been extensively studied [1]. The condition for 
CP symmetry violation in the SM requires at least thi^ee generations of fermions [2]. But, there 
is no upper limit on the number of generations from the first principle. In the SM the amount of 
CP violation is not large enough to explain the baiyon-antibaryon asymmetry in the Universe. The 
inclusion of fourth generation quarks leads to two extra CP phases in quark sector and possible 
larger CP violation [3], which is helpful for electroweak baryogenesis. With very massive quarks 
in the fourth generation, it has been proposed that the electroweak symmetry breaking may become 
a dynamical feature of the SM [4-7]. 

Heavy stable neutrinos with mass greater than ~ 1 GeV are possible candidates for the cold 
DM [8,9]. However, if the neutrino is the dominant component of the halo DM, the current DM 
direct search experiments have imposed strong constraints on its mass [10-15]. On the other hand, 
it is well-known that for a neutrino heavier than ~ niz/2, the cross section for its annihilation is in 
general too large to reproduce the observed DM rehc density. For the neutrino heavier than miy, the 
contribution from // channels decrease rapidly, but other channels such as Z^hP etc. are 

opened. For these processes the corresponding cross section does not decrease with the increasing 
of the neutrino mass, resulting in a relic density always decreases with the growing of the neutrino 
mass, and a thermal relic density far below the observed total DM relic density [16]. Thus the 
neutrino DM can only contribute to a small fraction of the relic density of DM and a small fraction 
of the halo DM density as well. Despite its very low number density in the halo, it can still be 
probed by the underground DM direct detection experiments due to its relatively strong coupling to 
the target nuclei, which provides a way to search for new physics beyond the SM complementary 



In this talk, we discuss the consequence of this possibility in a model with a fourth generation 
Majorana neutrino DM. The stabiUty of the fourth Majorana neutrino protected by an additional 
generation-dependent U{V) gauge symmetry which is anomaly-free. The details of the analysis can 
be found in Ref. [17]. We consider a simple extension of the SM with a sequential fourth generation 
and an additional U{\)f gauge symmetry. The U{\) extensions to the SM are well motivated from 
the point view of grand unification such as the 50(10) and E(, and have rich phenomenology [18] 
which can be reached by the on going LHC experiments. The flavor contents in the model are given 



All the fermions in the model are vector-like under the extra gauge interactions associated with 
U{\)f. The U{\)f chai^ges of the fermions could be generation-dependent. In order to evade 
the stringent constraints from the tree-level flavor changing neutral currents (FCNCs), the U{\)f 
charges Qqi for the first three generation quarks ai^e set to be the same, i.e. Qqt = Qq, (/= 1,2,3) 
while Qq4 = —3Qq for the fourth generation quarks. Similarly, the ^7(1)^ charges for the first 
three generation and the fourth generation leptons are Ql and —3Ql, respectively. In general, Qq 
and Ql can be different. For simplicity, we consider Qq = Ql = ^- With this set of flavor con- 
tents and U{l)f chai^ge assignments, it is straight forward to see that the new gauge interactions 
are anomaly-free. Since the gauge interaction of U{l)f is vector-like, the triangle anomalies of 
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[U{\)Ff, [SU{3)cfU{\)F and [gravity] ^[/(l)^ ai'e all vanishing. The anomaly oiU {\)y[U {\)f]^ 
is zero because the U{\)y hypercharges cancel for quarks and leptons separately in each genera- 
tion, namely + V) = ^nd K-^^l + J'w) = 0. The anomaly of [SU {2)lYU {1)f is also 
zero due to the relation Yfi= \ Qqi = and ^^^j Qu = 0. Thus in this model, the gauge anomalies 
generated by the first three generation fermions ai^e canceled by that of the fourth generation one, 
which also gives a motivation for the inclusion of the fourth generation. 

The gauge symmetry ?7 ( 1 is to be spontaneously broken by the Higgs mechanism. For this 
purpose we introduce two SM singlet scalar fields ^a,h which carry the U{\)f charges Qa = —2Ql 
and Qh = 62/, respectively. The U{X)f charges of d^ah are arranged such that 0^ can have Majorana 
type of Yukawa couplings to the right-handed neutrinos of the first three generations V/k (/ = 1,2,3) 
while 0/, only couples to the fourth generation neutrino Va,r. After the spontaneous symmetry 
breaking, the two scalar fields obtain vacuum expectation values (VEVs) ((^>a,i,) = Va.bj^- The 
relevant interactions in the model are given by 

7d 



— ^ij^ihUdiR — Y"jqiiHuiR — YfjijiHeiR — YijluHViR 
-\^u'WR^aVjR (/,7 = l,2,3)-^yf^(^),V4«-V((/),„(/)&,//)+H.c. (2) 

where /,• stand for left- and right-handed fennions, and H is the SM Higgs doublet. D^/,- = 
{d^ — igiZ"W^ — iYg2B^ — iQfgFZ'^)fi is the covaiiant derivative with the extra gauge boson 
associated with the U{1)f gauge symmetry, and gF the con^esponding gauge coupling constant. 
Since ^a,b are SM singlets, they do not play any role in the electroweak symmetry breaking. Thus 
Z' obtains mass only from the VEVs of the scalai^s 

'nl:=gUQlvl + Qlvl). (3) 

From the U{1)f charge assignments in the model, the four by four Yukawa coupling matrix is 
constrained to be of the block diagonal form 3 05 1 in the generation space. Since the U{1)f chai^ges 
are the same for the fermions in the first three generation, there is no tree level FCNC induced by 
the Z'-exchange in the physical basis after diagonalization. Thus a number of constraints from the 
low energy flavor physics such as the neutral meson mixings and the b ^ sy can be avoided. 

The direct search for the process e^e^ — >• Z' — )• t^l'' at the LEP-II leads to a lower bound on 
the ratio of the mass to the coupling to leptons: Mz' /gF > 6 TeV [19] for vector-like interactions. 



which corresponds to a more stringent lower bound: y Q^v^ + Q\v^ > 6 TeV. The cuixent seaixhes 
for narrow resonances in the Drell-Yan process pp ^ Z' ^ i^i^ at the LHC impose an alternative 
bound on the mass and the couplings of the Z' boson. For a model with sequential neutral gauge 
boson Zj^^ which by definition has the same couplings as that for the SM ZP boson [18], the latest 
lower bounds on its mass Mj, is 1.94 TeV from CMS [20] and 1.83 TeV from ATLAS [21] 
respectively. The bound on Mz'^^^ can be translated into the bound on the mass and couplings of 
the Z' in this model. For instance, gF < 0.029 for Mz' = 1 .44 TeV and < 0.005 1 for Mz' =0.94 
TeV, respectively. 

The fourth generation neutrinos obtain both Dirac and Majorana mass terms through the vac- 
uum expectation values (VEVs) of H and (ph- In the basis of (vl, v^)-'^ the mass matrix for the 
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fourth neutrino is given by 



mv = , (4) 

\ mo iiiM 



where mo = Y^vh/V2 with vh = 246 GeV and mm = Y^v^^j^fl. The left-handed components 
( v|™^ , v^^ ) of the two mass eigenstates are related to the ones in the flavor eigenstates by a rotation 
angle Q 

vIl' = -i{ceVL-S0V^), v^^' =seVL + cev'K, (5) 

where sg = sin 6 and cq = cos 6. The value of 6 is defined in the range (0, 7r/4) and is determined 
by 

tan20 = (6) 

niM 

with 6 =0 (7r/4) corresponding to the limit of minimal (maximal) mixing. The phase / is intro- 
duced to render the two mass eigenvalues real and positive. The two Majorana mass eigenstates 
are Xi = + v['^^'^ and Xi = ^21^ + ^2^'^, respectively. The masses of the two neutrinos are 
given by mi 2 = (y'm^H-4m^=FmM)/2. In terms of the mixing angle 6 they can be rewritten as 
m\ = {s0/c0)mD and m2 = {cg /se)mo- with mi < m2. Note that for all the possible values of 6 
the lighter neutrino mass eigenstate Xi consists of more left-handed neutrino than the right-handed 
one, which means that Xi always has sizable coupling to the SM Z" boson. Therefore the LEP-II 
bound on the mass of stable neutrino is always valid for Xi^ which is insensitive to the mixing 
angle. 

As the fermions in the first three generations and the fourth generation have different U{l)f 
charges, the fourth generation fermions cannot mix with the ones in the first three generations 
thr^ough Yukawa interactions. After the spontaneous breaking down ofU{l)f, there exists a resid- 
ual Z2 symmetry for the fourth generation fermions which protect the fourth neutrino Xi to be a 
stable particle if it is lighter than the fourth generation charged lepton 64, which makes it a possible 
dark matter candidate. 

The thermal rehc density of Xi is related to its annihilation cross section at freeze out. We 
numerically calculate the cross sections for ;^i;^i annihilation into all the relevant final states using 
CalHEP 2.4 [22]. In Fig. [l|, we show the quantity = Q.jf^jQ.DM the ratio of the relic density of 
X\ to the observed total DM relic density Q^oMh^ = 0.1 10 ±0.006 [23] as function of the mass of 
Xi for different values of the mixing angle d. The results show a significant dependence on the 
mixing angle 6. For smaller mixing angle 6 the couplings between X\ and gauge bosons W*,Z 
are stronger, resulting in a smaller relic density. The results also clearly show that due to the large 
annihilation cross section, Xi cannot make up the whole DM in the Universe. X\ can contribute to 
~ 20 — 40% of the total DM relic density when its mass is around 80 GeV. But for mi > m,, it can 
contribute only a few percent or less to the whole DM. 

However, since X\ has strong couplings to /i° and Z°, even in the case that the number density 
of Xi is very low in the DM halo, it is still possible that it can be detected by its elastic scattering 
off nucleus in direct detection experiments. Given the difficulties in detecting such a neutral and 
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stable paiticle at the LHC, there is a possibihty that the stable fourth generation neutrino could be 
first seen at the DM direct detection experiments. 

The generic formula for the differential event rate of DM-nucleus scattering per nucleus mass 
is given by 

dN _ PdmO,_^,^^^^ r-d\M, (7) 



where Er is the recoil energy, is the scattering cross section corresponding to the zero momen- 
tum transfer, mdm is the mass of the DM particle, = niDM^^N / {^nuM + "iw) is the DM-nucleus 
reduced mass, F{Er) is the form factor, and /(v) is the velocity distribution function of the halo 
DM. The local DM density Pdm is often set to be equal to po ~ 0.3 GeV/cm^ which is the local 
DM density inferred from astrophysics based on a smooth halo profile. Since the neutrino DM can 
only contribute to a small fraction of the relic density of DM, it is likely that it also contributes to 
a small fraction of the halo DM density, namely, its local density pi is much smaller than po. We 
assume that pi is proportional to the relic density of X\ in the Universe, namely 



Pi ^ ^Xi 



(8) 



DM 



or Tp « roi- Consequently, the expected event rates of the DM-nucleus elastic scattering will be 
scaled down by Vp. In order to directly compare the theoretical predictions with the reported ex- 
perimental upper limits which are often obtained under the assumption that the local DM particle 
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density is po, we sliall calculate the rescaled elastic scattering cross section 

d = rpO ^r^o, (9) 

which corresponds to the event rate to be seen at the direct detection experiments. Note that a 
depends on the mass of X\ through the ratio rp even when a is mass-independent. The spin- 
independent DM-nucleon elastic scattering cross section in the limit of zero momentum transfer is 
given by [24] 

where Z and A — Z are the number of protons and neutrons within the target nucleus, respectively. 
IXn = mim„/{mi +m„) is the DM-nucleon reduced mass. The couplings between DM and the 
proton (neutron) read 

fp(n)- L It, + (11^ 

q=u,d,s V q=c,b,t H 

with fj^^^ the DM coupling to light quarks and fp-^ = 1 — Y.q=ii4,s fj'q' ■ The coefficient in the 
model is given by 

2 m\mn 

ag=CQ^-^. (12) 

The value of aq is proportional to mi, thus larger elastic scattering cross section is expected for 
heavier Xi ■ The quark mass in the expression of cancels the one in the expression of . 
Thus there is no quark mass dependence in the calculations. 

In Fig. ^ we give the predicted spin-independent effective cross sections df^ for the fourth 
generation neutrino elastic scattering off nucleon as function of its mass for different values of 
the mixing angle 6. One sees that even after the inclusion of the rescaling factor rp, the current 
XenonlOO data can still mle out a stable fourth generation neutrino in the mass range 55 GeV < 
mi < 175 GeV which corresponds to r^^ < 1%. Thus the stable fourth generation neutrino must be 
heavier than the top quark, and can only contribute to a small fraction of the total DM relic density 
On the other hand, for m^, > 200 GeV, the cross section does not decrease with m^^ increasing, and 
is nearly a constant a,f' 1.5 x lO^'^'^cm^ in the range 200 GeV < m^^ < 400 GeV. This is due to 
the enhanced Yukawa coupling between the fourth generation neutrino and the Higgs boson which 
is proportional to m^^, as it is shown in the expression of a^. Similar conclusions are expected for 
other models in which DM particles interact with SM particles through Higgs portal, for instance, 
the singlet scalar DM in extensions of left-right symmetry model [25-28]. One can see from the 
Fig. ^ that the result is not sensitive to the mixing angle 6 either, which is due to the compensation 
of the similar dependencies on 6 in the relic density. For instance, the cross sections for the 
and Z^Z^ channel of XZ annihilation are proportional to Cg, which compensates the -dependence 
in the aq for the elastic scattering processes. 

The Majorana neutrino DM contributes also to spin-dependent elastic scattering cross section 
through axial-vector interaction induced by the exchange of the Z^ boson. At zero momentum 
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(GeV) 



Figure 2: Effective spin-independent cross section a^' which is a^^ rescaled by rp « ra for X\ elasti- 
cally scattering off nucleon as function of the mass of X\- Four curves correspond to the mixing angle 
= 10°(solid), 20°(dashed), 30°(dotted) and 40° (dot-dashed) respectively. The current upper limits from 
CDMS [29] and Xenon 100 [30] experiments are also shown. 



transfer, the spin-dependent cross section has the following form [24] 

afP = —Gj^f-^ {ap{Sp)+an{Sn))\ (13) 



where / is the spin of the nucleus, ap(„) is the DM effective coupling to proton (neutron) and 
(5'p(„)) the expectation value of the spin content of the nucleon within the nucleus. Gp is the Fermi 
constant. The coupling a^^^) can be written as 

-pin)=L^<"^ (14) 

where dq is the DIVI coupling to quark and A'y'"' is the fraction of the proton (neutron) spin carried 
by a given quai^k q. The coefficients in this model are given by 

du = -dd = -ds = ^. (15) 
V2 

For the axial-vector interactions, the coupling strengths do not depend on the electromagnetic 
chai^ges of the quarks. 

In Fig. |3| we show the predicted effective spin-dependent DlSl-neutron cross section d^^ as 
function of the neutrino mass for different mixing angles, together with various experimental upper 
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limits. Since a„ is independent of m■^^, the dependency of a^^^^ on the neutrino mass comes 
from the dependency of rp on nijf^^ , which can be seen by comparing Fig. ^ with Fig. |l|. The 
Xenon 10 data is able to exclude the neutrino DM in the mass range 60 GeV < m),^ 120 GeV, 
which is not as strong as that from the XenonlOO data on spin-independent elastic scattering cross 
section. For a heavy neutrino DM with mass in the range 200 GeV < m),^ < 400 GeV the predicted 
spin-dependent cross section is between 10^^° cm^ and 10^^^ cm^. 

In Fig. ^ we give the predicted spin-dependent DM-proton cross section Op^. The cross 
sections for Majorana neutrino DM scattering off proton and neutron are quite similar, which is 
due to the fact that the relative opposite signs in A„ and A^i are compensated by the opposite signs 
in di, and (i„. So far the most stringent limit on the DM-proton spin-dependent cross section is 
reported by the SIMPLE experiment [31]. The SIMPLE result is able to exclude the mass range 
50 GeV < m^^ < 150 GeV, which is compatible with the constraints from XenonlOO. In Fig. 0, we 
also show the upper limits from indirect searches using up-going muons which are related to the 
annihilation of stable fourth generation neutrinos captured in the Sun. The limit from the Super-K 
experiment is obtained with the assumption that 80% of the DM annihilation products are from bb, 
10% from cc and 10% from Tf respectively [32]. In the range 170 GeV < m^^ < 400 GeV, the 
limit from Super-K is ~ 5 x 10^^^ cm^. The IceCube sets a stronger limit dp^ < 2 x 10^'**' cm^ for 
the DM mass at 250 GeV [33]. This limit is obtained with the assumption that the DM annihilation 
products are dominated by W^W^. If the annihilation products are dominated by bb, the limit is 
much weaker, for instance a^f^ < 5 x 10"^^ cm^ for the DM mass at 500 GeV [33]. Note that in 
this model, the dominant final state is Z'^h^. The expected limit should be somewhere in between. 
Nevertheless, the IceCube has the potential to test these predictions. 

In summary, we have investigated the properties of stable fourth generation Majorana neutri- 
nos as dark matter particles. Although they contribute to a small fraction of the whole DM in the 
Universe, they can still be easily probed by the current direct detection experiments due to their 
relatively strong couplings to the SM particles. We have considered a fourth generation model with 
the stability of the fourth Majorana neutrino protected by an additional generation-dependent U{1) 
gauge symmetry. We have shown that the current XenonlOO data constrain the mass of the stable 
Majorana neutrino to be greater than the mass of the top quark. For a stable Majorana neutrino 
heavier than the top quark, the effective spin-independent cross section for the elastic scattering off 
nucleon is found to be insensitive to the neutrino mass and is predicted to be around 10^^^ cm^, 
which can be reached by the direct DM search experiments in the near- future. The predicted ef- 
fective spin-dependent cross section for the heavy neutrino scattering off proton is in the range 
10^'**' cm^ ~ 10^^^ cm^, which can be tested by the ongoing DM indirect search experiments such 
as IceCube. This work is supported in part by the National Basic Research Program of China (973 
Program) under Grants No. 2010CB833000; the National Nature Science Foundation of China 
(NSFC) under Grants No. 10975170, No. 10821504 and No. 10905084; and the Project of Knowl- 
edge Innovation Program (PKIP) of the Chinese Academy of Science. 
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